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' Abstract 
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■ We compute the complete two-loop 0{a) Wilson coefficients relevant for radiative 

decays of the B meson in the SM. This is a necessary step in the calculation of the 
0(aa" In" nib/Myy) corrections and improves on our previous analysis of electroweak effects 
in i? — > Xg'-f. We describe in detail several interesting technical aspects of the calculation 
and include all dominant QED matrix elements. In our final result, we neglect only terms 
originated from the unknown 0{aas) evolution of the Wilson coefficients and some sup- 
pressed two-loop matrix elements. Due to the compensation among different effects, we 
find that non-trivial electroweak corrections decrease the branching ratio by about 3.6% 
for a light Higgs boson, very close to our previous result. The corresponding up-to-date 
SM prediction for the branching ratio with > 1.6 GeV is (3.61 ± 0.30) x 10~^. 



1 Introduction 



The Standard Model (SM) calculation of the branching ratio for the inclusive radiative 
decay B Xs^ — BR^ in the following — has reached a high degree of sophistication 
(see ^ for a complete list of references and an up-to-date analysis) . Besides Leading Log- 
arithmic 0(a"L") (L = Innib/Mw) and Next-to-Leading Logarithmic (NLO) 0{a"L"^^) 
QCD corrections and non-perturbative Heavy Quark Effective Theory contributions, elec- 
troweak effects are known to play a non-negligible role In a previous work 0, we 

have considered in detail the electroweak corrections to this process, devoting special atten- 
tion to the interplay between QCD and electroweak effects. Photonic interactions generate 
logarithmically enhanced contributions which are suppressed by a factor a /as with respect 
to the QCD ones. The leading QED effects are therefore 0{aa^~^L^) and are known 
completely 0, while genuine electroweak corrections involving and W bosons start at 
the next order in the resummed logarithmic expansion. Ideally, one would like to have all 
these 0{aa'^L"') corrections under control. Since the 0{a) contributions to the coefficients 
of the four quark operators are all known 0,0, this would entail the following steps |@,|[: 

(i) the calculation of the two-loop 0{a) matching conditions for the magnetic operators 
Qj and Ql at a scale 0(Mvi,-); 

(ii) the QED-QCD evolution of the Wilson coefficients down to the B mass scale, includ- 
ing the calculation of the two and three-loop 0{aas) anomalous dimension matrix; 

(iii) the calculation of the one-loop and two-loop QED matrix elements of the various 
operators as well as of some yet unknown two-loop QCD matrix elements. 

Our analysis in 0| was based on the simplifying assumption that terms vanishing as 
Syir = sin 9yi.- can be neglected, unless they are enhanced by powers of the top mass Mt. 
In this case, introducing the SU(2)l coupling g and = 9^/4:71, all electroweak corrections 
are in fact O^awa^L^) or 0{aM^ /M^a^L^) and are included by step (i) only. Moreover, 
the calculation of the Wilson coefficients simplifies considerably. Although reasonable, this 
assumption should be verified — keep in mind that ~ 0.23. In particular, Z° boson 
corrections to the one-loop 6 ^ 57 magnetic penguin diagrams give rise to O(s^) terms 
which are not formally suppressed by an electric charge factor Q"^ = 1/9 or Qu\Qd\ = 2/9, 
unlike the purely QED corrections of steps (ii) and (iii). This happens, for instance, 
because of the mass difference between Z° and W bosons. Such O(s^) terms originate at 
the electroweak scale and affect only step (i). 

In this note we extend our calculation and compute the full 0{a) contribution to the 
Wilson coefficients of the b ^ s magnetic operators, thus completing step (i). The main 
difference (and technical hurdle) with respect to p[ is due to the presence of virtual photons 
in the two-loop SM diagrams. The resulting infrared (IR) divergences are removed in the 
matching with the effective low-energy theory of quarks, photons and gluons. Several 



2 



subtleties arise in the calculation, mostly linked to the presence of unphysical operators. 
This is explained in detail in Section 2, while Section 3 deals with the QED-QCD evolution 
of the coefficients and illustrates how 0{aa^L^) effects should be taken into account in 
the calculation of BR^. We also include all dominant 0{a) matrix elements and conclude 
reconsidering the SM prediction of BRy. 



2 The 0{a) matching 



Let us briefly recall the formalism. We work in the framework of an effective low-energy 
theory with five active quarks, photons and gluons, obtained by integrating out heavy 
degrees of freedom characterized by a mass scale M > Mw In the leading order of the 
operator product expansion the effective off-shell Hamiltonian relevant for the 6 — > 37 and 
b ^ sg transition at a scale /i is given by 
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Here Vij are the CKM matrix elements and Ci{fi), C^{^) and denote the Wilson 

coefficients of the following set of gauge invariant operators [P, PHTll 
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where e (qs) is the electromagnetic (strong) coupling constant, g^/j are the chiral quark 
fields, Ffj,i, (G^j,) is the electromagnetic (gluonic) field strength tensor, is the covariant 
derivative of the gauge group SU(3)c x U(1)q and T°- are the colour matrices, normalized 
so that Tt{T°'T'') = 6°'^/2. The s-quark mass is neglected in Eq. (|^) and in the following. 
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Notice that at the order we are going to work it is not necessary to consider the analogues 
of Qi and Q2 involving the w-quark instead of the c-quark. 

The above set of operators closes off-shell under QCD and QED renormalization, up 
to non-physical (evanescent) operators that vanish in four dimensions It consists of 

the current-current operators Q1-Q2, the QCD penguin operators Qs-Qe, the electroweak 
penguin operators Q7-Q10 and the magnetic moment type operators and Qf. It is the 
QED renormalization that forces us to introduce the operators Q7-Q10, in which the sum of 
the quark flavors is weighted by the electric charges Qg. The remaining six operators Qn- 
QiQ, characteristic of the process b ^ S'j (g), were chosen in such a way that they vanish 
on-shell up to total derivatives. Only operators of dimension five or six are retained. Higher 
dimension operators are suppressed by at least one power of m^/M^, while those of lower 
dimensionality can be removed by choosing suitable renormalization conditions in the full 
theory [|T2]|. In the present case this is achieved by requiring that all flavor off-diagonal 



quark two-point functions which appear at the one-loop level in the full theory vanish 
when the equations of motion (EOM) are applied, i.e. by using LSZ on-shell conditions 
on the external quark lines. We also use these renormalization conditions for internal 
virtual quarks in the full SM calculation and implement in this way a gauge invariant 0{a) 



definition of the CKM matrix 13 



In order to obtain the Wilson coefficients of the magnetic operator Q^, we calculate the 
off-shell amplitude 6 — > 57 in the full SM and in the effective theory at 0{a) and match 
the two results. Retaining only the leading terms in the off-shell amplitude in the 

full theory can be written in the following formQ 

= -^VtlVtb m bf^ , (3) 

where {s'y \Qi\b)^'^^ are the tree-level matrix elements of the operators in Eq. The 
perturbative expansion of the coefficients Ai reads 



A, = Af" + . (4) 



We calculate analytically the relevant one and two-loop amplitudes starting from the dia- 
grams generated by FeynArts 2.2 and retaining only terms which project on Qj after 
use of the EOM. All ultraviolet (UV) divergences in Aj^^^ are removed by electroweak renor- 
malization. We follow closely the procedure outlined in [@,|15|. The only additional ingre- 



dients not explicitly given in those papers are the right-handed down quark wave function 



renormalization (c^- = 1 — s^) 
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^In Eqs. (g|) and (|8|), the sum runs over Qi-Qw, Q7, Qf and, as we will explain later on, some evanescent 
operators. 
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and the complete 6-quark on-shell mass counterterm [xt = / M, 
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which had not been given in [0, . Notice that Eqs. (^) and (|^) do not include the 
photon contribution, for a reason that will become clear in a moment. The top mass Mt 
is renormalized on-shell as far as electroweak effects are concerned, while we use an MS 
definition at a scale /i for the QCD renormalization. 

We work in the background field gauge (BFG) This reduces the number of dia- 
grams to be considered. Moreover, if the electric charge is normalized at = q, as it is 
natural to do 0], its counterterm cancels identically against the background photon wave 
function renormalization factor, due to the BFG Ward identity |T^. The same holds in 
the case of an external gluon in the MS scheme. The regularization problems related to 
the definition of 75 in n = 4 — 2e dimensions are avoided as described in p| and we employ 
the naive dimensional regularization scheme with anticommuting 75 (NDR) throughout 
the paper. 

For what concerns the regularization of the IR divergences, we have adopted two differ- 
ent methods and found identical results for the Wilson coefficients. In the first method the 
IR divergences are regulated by quark masses (see |10|), while the second method consists 



in using dimensional regularization for both UV and IR divergences ||TT 



A second step involves the calculation of the off-shell amplitude in the QED effective 
theory. In general, we need effective vertices with both background and quantum photons. 
Interestingly, the latter introduce some gauge variant operators at 0{a.). In fact, on the full 
theory side there are heavy particle subdiagrams (see Fig. ||) that are coupled to quantum 
photons and contribute to gauge variant operators not included in the operator basis of 
Eq. (H). This is due to the gauge coupling of quantum photons with W and Goldstone 
bosons and is different from what happens in the case of the off-shell 0{as) matching 
l^,|Tl|. Indeed, at 0{as) only quark-gluon couplings and trilinear quantum-quantum- 
background gluon couplings are relevant and no gauge variant operator is induced. The 
appearance of gauge variant operators in the SM amplitudes is not surprising (see ||T9| 



for an example). 

We have explicitly verified that it is not necessary to take any gauge variant operator 
into account on the effective theory side. This follows from well-known theorems on the 
renormalization of gauge invariant operator^ [ITSj ^: gauge variant operators that mix 
with gauge invariant operators can be chosen so that they are all BRST-exact, i.e. they 

^The theorems apply to Yang-Mills theories, but extend to the full SM after imposing the anti-ghost 
equation |20|. 
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Figure 1: Subdiagrams involving the coupling of quantum photons 7 to the W and its corre- 
sponding Goldstone boson (left hand side) contribute to gauge variant operators (right hand 
side), as explained in the text. 7 denotes a background photon. 

can be written as the BRST- variation of some other operators, modulo terms vanishing 
by the EOM. Therefore, while gauge invariant operators generally mix into gauge variant 
operators^, the opposite is not true. Since we are eventually interested in the matrix 
elements of physical operators only, we do not need to include gauge variant effective 
operators in our basis. Of course this holds only as long as the regularization respects the 
symmetries, like it is in our case. 

In a similar way and because of the same theorems, the operators that vanish by 
EOM in the basis (0) do not mix into the physical operators of the same basis, and the 
renormalization mixing matrix is block triangular. This property drastically simplifies the 
computation at hand, as we will see in a moment. In particular, the renormalization mixing 
matrix Z is such that Zij = when Qi is EOM-vanishing and Qj is a physical operator. 

We have seen that effective vertices involving quantum photons are induced in the 
calculation. Even though contributions to gauge variant operators turn out to be irrelevant, 
the distinction is important in the case that the calculation is performed using quark masses 
to regularize IR divergences. For example, it turns out that the gauge invariant part of the 
off-shell b ^ s'-f effective vertex depends on whether the external photon is quantum or 
background. This can be explained by noting that the operators involving only background 
fields are combinations of truly gauge invariant operators and of operators containing also 
quantum gauge fields. This follows, e.g. from a decomposition of the kind p = p + ieQq^, 
where we used a hat to denote covariant background derivative, for the quantum photon 

■^At 0{a) the operators in our basis do not actually mix into gauge variant operators. 
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field, and Qq for tlie electric charge. The operators containing quantum gauge bosons 
eventually decouple from the calculation, as they are not gauge invariant. Because of the 
above decomposition, the coefficients of the operators involving only background fields are 
related to the coefficients of the operators in (H), as can be seen using Slavnov-Taylor 
identities of the kind used in p2 |. 

The effective theory calculation depends crucially on the IR regularization. The first 
method mentioned above (quark masses as regulator) can be applied at the diagrammatic 
level fl^. The effective theory diagrams are obtained by replacing hard (heavy mass) 



subdiagrams in the two-loop SM amplitudes with their Taylor expansions with respect 
to their external momenta. In principle, this method does not require a discussion of 
the effective operators. On the other hand, it is relatively complicated to implement. 
Here, we limit our discussion to the second method only, following In order to get 



the renormalized off-shell amplitude on the effective theory side, we need to reexpress 
Eq. (H]) in terms of renormalized quantities. The relations between the bare and the QED 
renormalized quantities are as follows 

Co = Z^e , mfe^o = ^mb"^fe , = Z^^A^ , go = 2"^^^ , Ci^ = ^ CjZji , (7) 



with Ze, Zm^, Zq and Zij the renormalization constant of the charge, the 6-quark mass, 
the quark fields and the Wilson coefficients, respectively. The relation between the renor- 
malization of the gauge field and of the electric charge is a direct consequence of the QED 
Ward identity. 

After renormalization the off-shell amplitude in the effective theory is given by 

Aff = -^^tlVtb J2 ^^^^^^ ^) ' (8) 



where Zj denotes a product of Z^, Z^ and Zq depending on the particular structure of the 
operator Qi and the Wilson coefficients may be expanded in powers of a as follows 



a(/x) = Cf)(/.) + f C«(/.) 



47r 



(9) 



As long as we are only interested in the Wilson coefficient of the magnetic photon 
penguin operator, it is sufficient to keep only terms proportional to (37 \Q'] \ h) in Eq. (0). 
Using the short hand notation (Qy) = (■57IQ7I&), the part of the off-shell amplitude in 
the effective theory needed for the matching of C7 is then written as 



Gf ^ 

71' 



Aff ^v:,ytb 



ZqZrrn, ^ CjZj^jj + Zq{Zm,, — 1) ^ CjZ_ 



i,13 



(10) 



where the second term proportional to Zj,i3 originates from the renormalization of the 
operator Qis. 
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Notice that the QED quark field renormahzation on the effective side can be avoided as 
it cancels in the matching against the photon contribution to the corresponding term in the 
SM. The same applies to the renormahzation of the 6-quark mass, which is retained only 
up to linear terms. Consequently, after checking the cancellation of the UV divergences in 
(|) we have omitted the photon contributions in Eqs. (||) and and simultaneously set 
Zrrif, and Zq to unity in the effective theory. This simplifies the following considerations. 

Adopting the MS scheme for the operator renormahzation the corresponding renormal- 
ization constants can be written as 

% = + £74' + ■ (11) 

The renormahzation constants Z^^^^ are found by calculating the UV divergent parts of 
Feynman diagrams in the effective theory. Within the scope of this computation, it is 



essential to carefully distinguish UV from IR singularities. As explained in Ref. [^, this 
can be done most easily by introducing a common mass parameter into all the propagator 
denominators including the photon ones. All renormahzation constants in the effective 
theory up to two loops are known from previous anomalous dimension calculations 
As we shall see later on, only five entries of the anomalous dimension matrix are relevant 
in the present computation and we have recalculated these elements to check the results 
mentioned above. Our results are in full agreement with Refs. |^,|^ and we will give the 
numerical values of the required renormahzation constants below Eqs. ( p!4D and (0). 

The last term in Eq. ([TT|) implies a finite renormahzation of Qi at zeroth order. Indeed, 
in situations where evanescent operators are present, the standard practice is to extend 
the MS scheme and to allow for a finite operator renormahzation. The finite terms zf^"* 
differ from zero when Qi is an evanescent operator and Qj is not, and their values are fixed 
by requiring that renormalized matrix elements of evanescent operators vanish in n = 4 
dimensions [p4|,p5|. This requirement also ensures that evanescent operators do not mix 
into physical ones Furthermore, in the case of the 6 — 57 calculation, it is well- 



known that some four quark operators can mix into the magnetic operators through 
one-loop diagrams at zeroth order in a and a^. Thus, not only we have finite terms in 
Eq. (pTD, but they appear at the lowest order in the coupling constant. 

The computation of the necessary matrix elements on the effective side is trivial, as 
we can set all the light particles masses to zeroQ. Accordingly, all loop diagrams on the 
effective side vanish in dimensional regularization, because of the cancellation between UV 
and IR divergences. Therefore only the tree-level matrix elements {Qi)^^^ are different 
from zero and higher order matrix elements do not play any role in the matching. Notice 
that due to the cancellation of UV and IR singularities the UV counterterms present in the 
tree-level matrix elements reproduce precisely the IR divergences in the effective theory. 
Furthermore, the IR divergence on the effective side has to be equal to the IR singularity 

We include only terms that are linear in the 6-quark mass. They originate from use of the EOM only. 
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on the SM side, to guarantee that the final results of the Wilson coefficients are free of IR 
poles. Eventually, all 1/e poles cancel out in C7, if the full and the effective theory are 
matched in the correct way. 

Bearing all this in mind, we are now able to extract from Eq. (|l^) those terms which are 
actually needed to calculate the 0(a) correction to the Wilson coefficient of the magnetic 
operator. First of all, we have to perform the tree-level matching by computing the relevant 
diagrams for the various operator insertions. Only C2, C7, C| and Cn-Cig are found to 
be non-vanishing at leading order. However, due to the triangularity of the mixing matrix 
the coefficients Cu-Ciq do not contribute to the first term in Eq. (|T0|) and therefore will 
not affect the matching conditions at the next order. Furthermore, as we set Zm^ equal to 
one also the term proportional to {Zm^ — 1) in Eq. ( |T0| ) does not contribute to at 0(a). 
Using Zq = 1 we thus obtain 
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(12) 



It is quite remarkable that, with the exception of the last term, only physical operators 
play a role in this expression, even though the calculation has been performed off-shell. 

The matching procedure between the full and the effective theory establishes the initial 
conditions for the Wilson coefficients a scale = 0{M^). Comparing Eqs. (D, (D and 
(|I^), the matching condition Afuiiifiw) = ^csi^J'w) translates into the following identities 
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(14) 



from which the Wilson coefficient of the magnetic operator up to 0(a) can be calculated. 
The leading order initial condition for the Wilson coefficient of Q2 is simply C2^\fiy^r) = 1 
and the elements of the mixing matrix needed for the next leading order matching of O7 are 



'2,77 



-58/243 and zi]]j^ 



8/9. Note that the renormalization constant ^2^^ ^^l^^^*^ to 



the mixing of the operators Q2 and Q^, is obtained from a two-loop calculation, as opposed 



to zi'^r^ which only involves a one-loop calculation. Whereas Z^^'j is regularization and 



(1) 



77,77 



renormalization scheme independent, ^2^7^ is scheme dependent. The value for 2^2 "77 given 
above corresponds to the NDR scheme — see 0. Notice also that in Eq. (14) the 0(e) 



77,77 



(1) 



terms of Cj^^^ yield a finite contribution when combined with the 1/e pole proportional 
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to Z^^-j . Indeed, the leading order matching needs to be performed up to 0(e). Exphcit 



formulas for the initial condition of C^*^'''' including 0(e) terms can be found in pTj 28|]. 

For what concerns the last term in Eq. (0), it is necessary to introduce the following 
evanescent operators 



[sLl^^bL) ^{ui^qL) + (1 + aie) - 
{sLl^fiL) J2^Q,{qLl^qL) + (1 + a2e) 



;i + a3e) 



(15) 



3Q7 - ^Q, 



where are arbitrary constants. In NDR inserting these operators into the one-loop 
b ^ sj penguin diagrams yields ZE^j-y = 4/9, ZE2,7-y = —4/27 and Ze^^^i-^ = 4/27. Notice 
that the last term in Eq. (0) does not depend on the special choice of evanescent operators 
adopted above, but it does depend on the choice of physical operators. For instance, in the 
operator basis of [0, all evanescent operators that project on Ql^ have vanishing Wilson 
coefficients, both at 0(a) and 0{as)- Therefore, in this basis evanescent operators do 
not affect the matching equations and is it not necessary to introduce them in Eq. (p^). 
Curiously, in the operators basis of Eq. (ET) the same holds only at 0{as). 



We have verified that all 1/e poles cancel in Eq. (14), and that the result for Cj ^ 



7(1) 
7,e 

coincides with the one obtained using quark masses for the IR regularization. In the 
latter case evanescent operators do not play any role in the matching, as their contribution 
to the matrix elements cancels against a corresponding term stemming from the finite 
renormalization z'^^\^. 

Let us now turn to the matching for the Wilson coefficient of the chromomagnetic 
penguin operator Qg. The calculation for the b s gluon off-shell amplitude proceeds in 
the same way as above. Adopting the notation (Q^) = {sg \Ql \ b), we see that the analogue 
of Eq. (PD is 
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from which we obtain 



9(0), 



(17) 



7(1) ^(0) 
^2,89'-"2 



7(1) r^9(0)/,, 
^89,89*^^8 lA^i- 



(0) ^E{1) 



(/i 



Wj 5 



10 



where Z. 



(1) 

2,8g 



-23/81, Z. 



(1) 

77,8g 



-4/3 and Z 
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4/9. The renormahzation constants 



which describe the mixing of evanescent operators into physical ones read ZE^^sg = ~4/3, 
ZE2,8g = 4/9 and Ze^^^sq = —4/9. Again, all IR poles cancel in Cf and the result coincides 
with the one obtained with the other method. 

We now recall that the relevant quantity entering the calculation of BR^ is not Cy (/Xb) 
with Hb = 0{mb) but a combination Cy'*^ (/ife) of this Wilson coefficient and of the coeffi- 
cients of the four quark operators. This combination is the coefficient of (Q7) '^'^'^ calculated 
on-shell. It follows from this definition that, unlike Cy, the effective coefficient is regu- 
larization scheme independent at LO |^ and does not depend on the basis of physical 
operators. In NDR the two combinations relevant for B —>■ Xgj and B XsQ are 
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The 0{a) contributions to the Wilson coefficients at /i = M„ in our operator basis Eq. 
can be found from those in the operator basis of after a basis transformation which in 
four dimensions is simply 



(20) 



where C are the Wilson coefficients in the basis of 0]. They are given in Eqs. (8.111)- 
(8.117) of that review. The matrix TZ is the extension of the same matrix of ||2^ and is 
needed only for the physical operators Qi-Qio, Qj and Qf: 
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In fact, beyond the leading order, the operator basis must be supplemented by a definition 
of the evanescent operators. This definition corresponds to a choice of scheme and it is 
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different, for instance, in the standard basis of 0] and in the operator basis of P,p9|,p0 
On the other hand, a change of scheme can be in general accommodated by an additional 
non-linear term in the transformation Eq. pO) 



C7(/i) 



1 + 



a 



(22) 



where Af^ and Afe are matrices that depend on the way the projection on the space 
of physical operators is implemented in the effective theory calculation (which in turn 
corresponds to a definition of evanescent operators). At the order we are interested in, 
they affect only Ci and C2. 

In the following, for definiteness, we follow the convention of |2^,|3^, whose basis of 
physical operators is a subset of Eq. (|). Recalling that ci]l{My^) was obtained in [Q in the 
standard basis, we have calculated the matrix Afe that connects the two different schemes. 
As a result, the non-vanishing 0{a) contributions to the Wilson coefficients of the four 
quark operators and of the evanescent operators at yU = M^y are given by 
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Figure 2: Electroweak corrections to the Wilson coefficients (^^(M^) and C|(Mvv). The 
dashed lines represent the results of with their error estimates, the solid lines the complete 
corrections to the Wilson coefficients at My^. 

Here we have left explicit the extra scheme dependent term 1/9 in C2^l{M^): it is numer- 
ically very small. 

The final results for Cj^^^ and Cfg^"* are quite lengthy. We give instead two accurate 
approximate formulas for the O^g"^) contributions to C^'^^ {M^) and C|''^^(M,^), which are 
valid when the effective Hamiltonian is normalized in terms of G^? as in Eq. (P: 



g,eff{l) 



1.11 - 1.15 1 



Ml 
1702 



- 0.444 In- 



H „ , 2 



100 



0.21 In' 



100 



-0.513 n— n — - 
100 170 



104\ 



+ ( ^-Cf^ - ^ 1 In 




-0.129 In 



243; Ml ' 

— - 0.0244 In^ — 
100 100 



-0.037 In 



9 ^ 3 ^ 



8lJ Ml 



(25) 



Here M^ is the Higgs boson mass expressed, like Mf in GeV. In Eq. we use the coupling 
ce{fiw) ~ 1/128, while in general we employ = 0.23, corresponding to g'^ = 4:\/2GfMI, 
Myy = 80.45 GeV and = 91.1875 GeV. Eqs. (|2|) reproduce accurately_[within 1.5%) 
the analytic results in the ranges 80 GeV < M^ < 300 GeV and 160 GeV <Mt< 180 GeV. 
We stress that Eqs. (^) are independent of the choice of the scale fit in the QCD top mass 
definition: it is sufficient to calculate Mt{fit) and employ it in Eqs. (pS]). Different choices 
of fit lead to different NLO QCD corrections, but they are higher order effects as far as 
the present calculation is concerned. The iiyy dependence of the effective coefficients agrees 
with ill. 
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The size of the electroweak corrections to C^'*' and C|'*^ relative to the one-loop 
results is shown in Fig. as a function of the Higgs mass. To compare directly the results 
in Eq. (pSj) with the approximate ones in 0, we have used the same central value Mj = 
175.5 GeV in the plots. First, notice that the Higgs mass dependence is identical, as 
should be expected since all the diagrams involving the Higgs boson also involve a charged 
boson. Therefore, these diagrams are not sensitive to the Z-W mass difference or to 0(s|,) 
couplings. Numerically, we see from Fig. |^ that the difference is larger than estimated in |Q . 
Although an expansion of the results in powers of converges quickly, it turns out that 
its second term of O(s^) is larger than naively expected, and that the two-loop correction 
is very sensitive to the M^-Mi^ difference. 



3 QED-QCD evolution and the decay B Xg^ 



The relevant quantity in the evaluation of BR^ is the effective Wilson coefficient of the 
magnetic operator at a scale fib ~ ^b- In the resummation of QED and QCD logarithms 
one usually keeps only terms linear in the electromagnetic coupling a, whose running is also 
neglected. In this case, the general structure of the evolution at 0(aa"L") is well-known, 
see for instance Using the notation = (Ci, . . . , Cio, C]''^ , Cf'*^ ) and restricting to 
the physical (on-shell) operators in Eq. (|]), the coefficients at a scale are given in terms 
of the coefficients at the scale M^^ by 



where 

f/(/i, M^r, a) 



An L 



(26) 



(27) 



The first two terms give the pure QCD evolution. The matrices f/*^*-* and R^^^ are deter- 
mined by the anomalous dimension matrix of the operators in question and by the QCD 
P function. Explicit expressions for U^^\ U^^\ and R^^^ can be extracted from p3|. 
R^^^ is presently unknown: it requires the evaluation of the two and three-loop anoma- 
lous dimension matrix at 0{aas). From the point of view of the expansion in as in the 
renormalization group improved perturbation theory, U^^^ and U^-^^ are 0(1) and 0{as), 
respectively. R^^^ and R^^'^ are 0(1/ as) and 0(1), respectively. Expanding in as and a we 
obtain 



0(/i) = f/(°)(^,M, 



'^(0) ^ Qs(^w) ^(i) 

4:71 " 



+ 



a 
Atx 



f/(°)(/i, M^)0(^) + i?W(/i, M^) + '^AMAcm\ + M^)o(°) 



14 



All the coefficients on the right hand side are understood at the scale Mw The first line 
results from pure QCD evolution. The second line mixes QED, electroweak, and QCD 
effects. After the calculation of the two missing elements of (7e^^ in Section 2, the only 
unknown part of Eq. (^) relevant for (chromo) magnetic decays is the last term. 

Turning to the particular case of C^''^^ and neglecting the unknown term M'^^fi, Mvir)C^^\ 
we see from Eq. (|^) and that the 0{a) terms are given by 



(29) 



The first term on the right hand side corresponds to IJ^^\^, M^)Ce in Eq. ( p8D and takes 
the form 



(30) 



- (0.448 - 0.49 r^) C^]l{M^) + (0.362 - 0.454 77) C^'i(M^) 

+ (5.57 - 5.86 r/) Cg(M^) + (0.321 - 0.47 r^) Cg(M^) 

+ (1.588 - 2.89 r/) eg (M^), 

where rj = as{M^y) / as{fib) ~ 0.56 for /i^ = m^. Here we have given analytically only 
the cofactors of Cy'*'^ and C^'*^^ . The other terms are more involved and are given in an 
approximate form, valid within 1% for values of 77 between 0.5 and 0.6. However, they can 
all be easily determined from the anomalous dimension matrices given in The last 
four terms have been given in 0] in the operator basis of 0. 

The second term in Eq. (^ corresponds to R^^\fi, M^){C^°^ + as{M^)/{4n)cP) in 
Eq. (pHF) and is given by 



in 



«s(M6) 



88 
575' 



~ I 1449 ' 



704 
1725' 



-0.0449 + 0.2504r7 - 0.236//^ 
-(0.381 - 0.556 r/)r7Cf](M^, 



+ (0.15-0. 178 r7)r7C{^](M, 



(31) 



The 0{as) coefficients C^^]{Mw) can be found in |T^,^,^. The approximate expressions 
are valid within 1% for 0.5 < r] < 0.6. 

We are now ready to give a numerical value for the 0{a) Wilson coefficient at yUj, = 
4.7 GeV using Eq. (0). Renormalizing the top mass at = Mt = 165 GeV, we find 



C?f^'H/ife) =4.172 -1.312 In 



100 



0.615 In^^ + 2.360 



100 



(32) 



^ Table 2 in H allows to change from the basis of M to that of 
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where the first three terms correspond to Cj'jf and the last one to C^'^- Notice that for 
a hght Higgs boson the first term (formally 0(aa"L")) is twice the second one (formally 
0(aa"~^L")). We interpret this as evidence that purely electroweak 0{aw) effects are 
dominant with respect to purely QED effects. 

To see how electroweak corrections affect the calculation of BR^ it is sufficient to recall 
that, for fib = nib, the perturbative QCD expression for the b — > sjl^g) decay is proportional 
to 



2 

+ B{Eo), (33) 



where Eo is the maximal photon energy in the 6-quark frame and B{Eq) originates from 
bremsstrahlung diagrams. Tg and re originate from the 0{as) and 0{a) matrix elements 
of the physical operators. Tg has been computed in with the exception of ra^s, • • • , ''"io,s- 
It is easy to see from these papers that^ 

2 1 11 

while rgg e = 0. Numerically, the effect of r2,s and rj^^g in the calculation of the inclusive 
branching ratio is quite important. On the other hand, the 0{as) contribution to B{Eq) 
changes BR^ by less than 4% if 1 GeV < Eq < 2 GeV. We therefore conclude that the 
only potentially relevant QED matrix elements are the virtual corrections parameterized 
by Ti^ei i^2,e and rr^y^e- Hence, in our numerics we will neglect the unknown last term in 
Eq. (pSD, the unknown QCD matrix elements (which are in any case suppressed by small 
Wilson coefficients), and the remaining real and virtual QED contributions to the matrix 
elements. It has been observed in that splitting the charm and top quark contributions 
in Eq. (|33D and normalizing them in an asymmetric way leads to an improved perturbative 
QCD expansion. In the evaluation of the electroweak corrections, however, this would be 
an unnecessary complication. 

We stress that, as we neglect the last term in Eq. ([28| ) and some contributions to the 
matrix elements, our evaluation of 0{aa'^L"') effects in B Xs'~f is incomplete, although 
we are confident that it should provide a good approximation. Our numerical result is valid 



in the NDR scheme supplemented by the definition of evanescent operators of ||29|,pO|]. An 
analysis of the way the scheme dependent terms recombine can be found in . The scheme 
dependence of our result is introduced in Eq. (^) and in the rj^e- The one from Eq. (|23|) is 
numerically negligible (less than 0.01% on BR^). All the residual scheme dependent pieces 
would be cancelled by corresponding terms in the anomalous dimension matrix at 0{aas), 
if it were available. 

Additional 0(a) contributions are introduced by normalizing BR^ in terms of the 
semileptonic branching ratio, as there are well-known QED corrections to the semileptonic 



^Because of the definition of B{Eq) adopted in we use here rj^^s — |(4 — tt^) as given in |Q. This 
convention is different from the one of fSOfl . 
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a 
An L 



decay amplitude Unfortunately, only the leading logarithmic term is known. The 

final 0{a.) contribution to the expression under absolute value in Eq. (|33|) is therefore 

+ r,,,Cr^^'^ (m,) - 4 Cr^'^ (m,) In ^1 . ^^^^ 

Using the reference value M„ = 115 GeV and mc/mf, = 0.22 as in [jl|, we find numerically 

£ew = 0.0025 + 0.0014 + 0.0004 + 0.0028 = 0.0071 , (36) 

which updates Eq. (4.6) of |[l|]. Here the first and second terms correspond to the U and 
R components of C^'^^^^\fib) (in H] they were 0.0035 and 0.0012, respectively). The third 
term derives from the QED matrix elements and was not included in the paper mentioned 
above. The last term, 0.0028, is due to the QED corrections to the semileptonic decay am- 
plitude and is the same as in Notice that the first term, although formally suppressed 
with respect to the second one, is larger, as it incorporates all purely electroweak contri- 
butions. The total effect of the QED and electroweak corrections in ^ew on the branching 
ratio is a 3.6% reduction while the 0(Q;a"L") contributions alone lead to a 1.6% reduc- 
tion. As different contributions accidentally compensate each other, ^ew is almost exactly 
the same that was used in |]l|. Incorporating all perturbative and non-perturbative QCD 
corrections and using the same numerical inputs as in [jl| , we therefore obtain for different 
values of the cutoff photon energy in the B meson frame the following results 

BR [B ^ X,7] = (3.74 ± 0.30) x 10"^ , 

BR [5 ^ X.7]^^>,.eGev = (3-61 ±0.30) X 10-\ 

which are very close to those given in the paper mentioned above. Here the errors are 
estimates of theoretical errors also based on the analysis of [0]. One can compare the first 
of these two results with the present experimental world average BR^ = (3.23±0.42) x 10 



-4 



4 Conclusions 

We have calculated the complete 0{a) Wilson coefficients relevant for radiative weak de- 
cays and described the implementation of 0{aa^L"') effects in detail, including also the 
dominant QED matrix elements. The final impact of these contributions on the branching 
ratio of -B — s> Xg'y is roughly 1.6% for a light Higgs mass ~ 100 GeV, and decreases 
slowly for larger values of M^. 

We have discussed in detail the role played by unphysical operators in the calculation. 
We have adopted two different methods to regulate the IR divergences and clarified the 
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subtleties that arise in the two cases. In contrast to the off-shell 0{as) calculation |TO| , p 
evanescent operators turn out to play a crucial role in the 0(a) computation. We have 
also explained the relevance of gauge variant operators in our calculation. 

Our results improve upon existent calculations and put electroweak corrections 

to B Xsj on a firmer basis, although numerically the change is negligible. The depen- 
dence of CY'^^{fib) on heavy degrees of freedom is now completely known at 0(a). We 
have also included the dominant 0(a) matrix elements. Still, not all the 0{aa^L") con- 
tributions to radiative decays are under control. The uncalculated corrections are related 
to the QED-QCD evolution (last term in Eq. (^)) and to some suppressed QED matrix 
elements. As the electroweak contributions to QCD and electroweak penguin operators 
are relevant to our discussion, some two-loop QCD matrix elements are also still missing. 
The incompleteness of our calculation makes it scheme-dependent, but, as we have noted 
above, the scheme dependence is remarkably small. On the other hand, the calculation 
of the missing contributions would require a significant effort. In the meanwhile, we note 
that: (i) the leading 0(aa"~^L"') corrections affect the branching ratio of 5 — > only in 
a minor way (—0.6%) and (ii) the QED matrix elements are very small, although formally 
of the same order of the matching corrections. Therefore, one might expect the missing 
subleading QED effects to be eventually small. 



We are grateful to Antonio Grassi and Mikolaj Misiak for many helpful discussions and 
to Andrzej Buras for a careful reading of the manuscript. 
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